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I. INTRODUCTION 



The quantum mechanical systems could be investigated in the view of two basic points. 
One of them is the studying of bound states to handle the necessary information about the 
system under consideration. The other point is solving of scattering problem for a given 
quantum mechanical system under the effect of a potential. So, one has to study both of 
bound states and scattering states of a quantum mechanical system under consideration to 
achieve a complete information about it. Some efforts have been made about the scattering 
problem for a relativistic and/or non-relativistic system under the influence of different types 
of potentials, such as Manning-Rosen potential [1, 2], Eckart potential [3, 4], Poschl- Teller 
potential [5], Hulthen potential [6], Woods-Saxon potential [7-9], cusp potential [10], and 
Coulomb potential [11]. The scattering problem in the case where the mass depends on spa- 
tially coordinate has become a particular part of that problem, and has been received great 
attention to study of scattering states for a given quantum system [12-15]. The position- 
dependent mass formalism is a useful ground to explain the electronic properties of quantum 
wells and quantum dots [16], semiconductor heterostructures [17], and impurities in crystals 
[18-20]. 

In this paper, we solve the following one-dimensional Schrodinger equation {h = 1) 

obtained from the Hamiltonian [14] 

H^Up-p)+V. (2) 
2 \ m / 

for the Woods-Saxon potential to study the scattering states within the framework of 
position-dependent mass formalism. The effective-mass Schrodinger equation could be trans- 
formed into Heun's equation [21] which is a Fuchsian-type equation with four singularities 
[14] by using a coordinate transformation. We obtain the wave function in terms of Heun's 
function and then we find transmission and refiection coefficients by studying the asymp- 
totic behavior of the wave function at infinity. We write also the transmission and reflection 
coefficients for the case of constant mass by using the properties of Heun's function and also 
the continuity conditions of the wave function at x — 0. We flnd the wave function for the 
case of constant mass in terms of hypergeometric functions and plot the wave functions for 
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completeness. In nuclear physics, the Woods-Saxon potential is used to construct a shell 
model to describe the single-particle motion in a fusing system [22] and the potential plays 
an important role within the microscopic physics because of describing the interaction of a 
nucleon with a heavy nucleus [23]. 

The work is organized as follows. In Section II we obtain exactly scattering state solutions 
of the Woods-Saxon potential and transmission and reflection coefficients in the case of 
position-dependent mass. We study also the same quantities in the case of constant mass. 
The conclusions are given in Section III. In Appendix A we list some equalities related with 
Heun's function required for this work. 

II. SCATTERING STATE SOLUTIONS 

The Woods-Saxon potential has the form 

and we parameterize the mass function as 

m{x) = (mo - mi) (M - y^^) , (4) 

where M — (mo -l- mi) /{mo — mi) and Vq, S, mo and mi are positive parameters. The form 
of the mass function is strongly similar to that of the potential. We could exactly solve the 
problem because of this form and also study the results for the case of constant mass. By 
using the transformation y = (1 + e^^)~^ and inserting Eq.(4) and Eq.(3) into Eq.(l), we 
obtain the differential equation (0 < y < 1) 

\y y-i y-Mj 

where 

al = {2/S^){mo-mi)Vo; -aj = {2/6^){mo + mi)E ; -aj = {4/6^)mi{E + Vo) . (6) 
To obtain a Fuchsian-type differential equation from Eq. (5), we use a new transformation 

i^{y) = y^Hy-'^rfiy), (7) 
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which gives a Heun's-type equation given as in Eq. (Al) in Appendix A 

\ y y-i y-Mj 



(8) 

The general solution of Eq. (8), which is regular in the neighborhood of y = 0, is written in 
terms of the Heun's function as [14] 

f{y) = AH{M, -[-aj + (a2 + as){l + 02 + a3)]M + 02; 02 + as - ai, 

02 + 03 + oi, 1 + 2a2, , (9) 

where the constant A will be determined below. 

Let us first investigate the limit x ^ oo[y e~^^ O), which gives /(O) = ^4 in Eq. (9) 
and the solution '^{y) Ay"-^ — Ae~^'^^^ becomes 

ij{x) = Ae-'^'^ , (10) 



where ki — ^/2(mo + mi)E and we have used the property of H{a, b; a, 7, S; 0) = 1. 

To study the behavior of the solution Eq. (9) for x — )■ —oo{y — )> 1), 1 — y ~ e^^, we use 
Eq. (A5) of Appendix A, which changes the argument y to 1 — y. Thus, we obtain the 
Heun's function in Eq. (9) as 

H{M, -[-al + (a2 + a3)(l + 02 + a3)]M + 02; 02 + 03 - ai, 02 + 03 + Oi, 1 + 2a2, = 
DiH{l - M, [-a? + (02 + a3)(l + 02 + a3)]M + a? - (02 + 03)^ - 02; 
02 + 03 - ai, 02 + 03 + ai, 1 + 203, -1; 1 - y) + £>2(1 - ?/)~^"' 
xif(l - M, [-ai + (02 - as){l + 02 - a3)]M + - (02 - 03)^ - 02; 

02 - 03 + oi, 02 - 03 - oi, 1 - 2a3, -1; 1 - y) , (11) 
where the constants Di and D2 are written by using Eq. (A6) in Appendix A 

Di = H{M,-[-(4 + {a2 + a3){l + a2 + a3)]M + a2; 

02 + 03 - «i, 02 + a3 + oi, 1 + 2a2, -1; 1) , (12) 
D2 = H{M, -[-al + (02 - a3)(l + 02 - a3)]M + 02; 

02 - 03 + ai, 02 - 03 - oi, 1 + 2a2, -1; 1) . (13) 



Using Eq. (11) we obtain the solution in Eq. (9) as 

V'(y) ^ >l(-l)"H^ie^"''' + 026-^"'^} , 



(14) 



which gives 



(15) 



where /c2 = \/A.mi{E + Vq) and we set ^4 = (—1) "^/-Di. Thus, we achieve the following 
form of the wave function for the limit x icxo 



(16) 



As a result, we recover the asymptotic behavior of a plane wave coming from the left-hand 
side. 

We can write the wave function explicitly 



V'(x) = (_l)2a3(l^gfe^-(a2+a3)g5a3X 

H{M, b + 02; 02 + 03 - oi, 02 + 03 + oi, 1 + 2a2, -1; 



X 



(17) 



H(M, 6 + 02; 02 + 03 - oi, 02 + 03 + oi, 1 + 2a2, -1; 1) 

where b = —[—al + (fl2 + «3)(l + 02 + 03)]-^- Finally, we give the reflection and transmission 
coefficients for the case of position-dependent mass, respectively 



\R\ 



H(M, b' + 02; a2 — as + oi, 02 - 03 - oi, 1 + 2a2, -1; 1) 



H{M, b + 02; a2 + as — ai, 02 03 + oi, 1 + 2a2, -1; 1) 
where b' — —[—a? -I- (02 — a3)(l + 02 — a3)]M, and 

1 



(18) 



ITP = — 



(19) 



k2 \H{M, b + a2] a2 + a^- ai, a2 + as + ai, 1 -h 2a2, -1; 1)|^ 

In order to investigate the dependence of the reflection coefficient to the energy E, we 
rewrite Eq. (18) in the following form by interchanging a •<-> ^ in Heun's function 

H{M, b' + 02; 02 - as ai, 02 - 0.3 - o-i, 1 2a2, -1; 1) 



\R? - 



X 



H{M, b + a2] a2 + a^. + ai, a2 + - ai, 1 + 2a2, -1; 1) 
H{M, b' -\- 02; 02 - 03 - oi, a2 - 03 ai, 1 2a2, -1; 1) 



H{M, b 02; 02 03 - oi, 02 03 oi, 1 2a2, -1; 1) 



(20) 
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By using Eq. (A7) in Appendix A and keeping in mind that aj — a\/M — a\/{M — 1), Eq. 
(20) gives 



\R\ 



X 



X 



[(M - l)a2 + Ma^f 
[(M - l)a2 - Ma^Y 
H(M, 6' + 03 - oi; 02 - 03 + oi, 02 



03-01 + 1,2 + 202,0;!) 



H(M, 6 - 03 - Oi; 02 + 03 + oi, 02 + 03 - oi + 1, 2 + 2o2, 0; 1) 
H{M, 6' + 03 + oi; 02 - 03 - oi, 02 - 03 + oi + 1, 2 + 2o2, 0; 1) 



21) 

i7(M, 6 - 03 + Oi; 02 + 03 - oi, 02 + 03 + oi + 1, 2 + 2o2, 0; 1) ' ^ 
This equation enables us to analyze the dependence of reflection coefficient to the energy E 
when the energy goes to infinity. In this case, Eq. (21) gives 



2_ / V2mi - ^niQ + nil 



H{M, -M{a2 



03)^; 02 



03,02 - 03,202,0; 1) 



if(M, -M(o2 + 03)2; 02 + 03, 02 + 03, 2o2, 0; 1) 



(22) 



Using the equality if(o, 6; a, /3, 7, 0; y) = 2-^1(0;, /3; 7; y) (for b — —aa/3) [14] and also 



2^1(0;, a; 7; 1) ^-a. 



,7— >CX) 



, we obtain 



\R\ 



E^oo I 



/ \/2mi — ^ynl^^'+n^ 



Eq. (23) shows that the reflection coefficient increases up to the value obtained in Eq. (23) 
while changing with energy. Fig. (1) shows the variation of the reflection and transmission 
coefficients as a function of the energy E in the position dependent mass case. In Fig. (2), 
the effect of the mass parameters mo and mi on the reflection and transmission coefficients 
are given. It is seen that the reflection coefficient decreases linearly with mass parameters 
while the transmission coefficient increases with the growing values of the parameters. In 
the Figs. (1) and (2), It could be seen that the unitarity condition |i?P+ |rp = 1 is satisfied 
in the constant and position dependent nass cases. In Fig. (2), we see that the reflection 
coefficient can not take zero value for the case of < Vq which is agreed with quantum 
mechanical results. 

Now, we begin to give the results for the case of constant mass, which means that mo = 
mi, starting from the wave function. With the help of Eq. (A8) in Appendix A, we write 
the wave function 



mo=mi 



X 



2^1(1 + 02 + O3, O2 + O3; 1 + 2O2; 



2^1(1 + O2 + O3, O2 + O3; 1 + 2O2; 1) 



(24) 
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which can be written in terms of Gamma functions 



mo=mi 



Sx\—(a2 



+a,)^5a3x '^{(^2 " a3)r(l + 02 - Q3) 



r(l + 2a2)r(-2a3) 



X 2-^1(1 + 02 + a3, a2 + 03; 1 + 2a2; 



1 + e 



5x 



)• 



(25) 



where used the relation of hypergeometric function 2-^1(0;,^; 7; 1) = Y^(J^-a)r(^-p) ■ "^^^ P^" 
rameters given in Eq. (6) in the case of constant mass become (mo = mi = m) 



-al^ {A/S^)mE ; -aj = {4:/S^)m{E + Vo) 



(26) 



We depict the wave function for two different values of parameter sets in Fig. (3). It is 
seen that the wave function exhibit an oscillatory behaviour for x < and exponentially 
decreasing in the region x > 0. The oscillating behaviour of the wave function given in Eq. 
(16) for a; < is a purely quantum mechanical interference effect between the incident and 
reflected waves [23]. The wave function in the region x > goes to zero due to the potential 
given in Eq. (3). 

We give the reflection and transmission coefficients for the case of constant mass. Using 
Eq. (A8) in Appendix A and the relation 2^1(0;, ^; 7; 1) = r(7-i)r(7-^) (1^)' ^® 

obtain 



\R\ 



mo=mi 



2Fi(a2 - a3 + 1, a2 - 03; 1 + 2a2; 1) 



2-^1(02 + 03 + 1, 02 + 03; 1 + 202; 1) 



r(2a3)r(a2 - a3)r(a2 -03 + 1) 



r(-2a3)r(a2 + 03)1(02 + 03 + 1) 



(27) 



and similarly from Eq. (19) 



in 



mo=mi 



ki 1 

h I 2-^1(02 + 03 + 1, «2 + 03; 1 + 2a2; 1) 



k2 



r(a2 - a3)r(a2 - ^3 + 1' 



r(l + 2a2)r(-2a3) 



(28) 



It should be noted that we must apply the continuity condition to obtain a relation 
between the coefficients written in Eq. (16). The condition that the wave function and its 
derivative must be continuous at x = gives A;2(l — |-Rp) = /ci|T'p [24, 25]. In Fig. (4), we 
plot the variation of the reflection and transmission coefficients according to the energy E in 
the case of constant mass. The reflection coefficient goes to zero when the energy increases 
while the transmission coefficient goes to unity. It could be interesting to study the limiting 
case of 5 — > 00. In that case the potential function becomes y(x) — > and the mass function 



goes to 2m. It means that the reflection and transmission can not appear (Eqs. (27) and 
(28)) as expected. In addition, in the hmiting case 5 — > — oo we obtain a step potential from 
Eq. (3) and Eq. (4) gives us m{x) — >■ 2m. Thus, we get the reflection coefficient as 



where ki — \/ AmE and k2 — ^/4m{E + Vq). 

III. CONCLUSION 

We have exactly solved the one-dimensional effective mass Schrodinger equation for the 
Woods-Saxon potential. We have found the wave functions in terms of Heun's function. The 
reflection and transmission coefficients are calculated by using the asymptotic behaviour of 
the wave function at infinity To analyze these coefficients in the case of position-dependent 
mass, we calculate the reflection coefficient in the limit E ^ oo. They are plotted as a 
function of mass parameters in Fig. (2). One can see that the unitarity condition in the 
scattering problem given asli^p-l-lTp^lis satisfied in the position dependent mass case 
also. We have also obtained the wave function, reflection and transmission coefficients in 
the constant mass case. They are presented in the Figs. (3) and (4). 
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Appendix A: Useful Equalities of Heun's Function 

Heun's equation with the following form 




(29) 



y-ajdy y{y - l){y - a) 



6 \ d a^y + h 



f{y)^0, (Al) 



has a solution in the neighborhood oi y = 



f{y) = H{a,b;a,p,^,S;y) , 



(A2) 



8 



and two linearly independent solutions in the neighborhood oi y — 1 [14] 

f{y)^H{l-a,-b-a/3;a,/3,l + a + l3-^-5,5;l-y), (A3) 

and 

f(y) = (1 - yy+^-^-PH{l - a, -6 - a/3 - (7 + (5 - a - /3)(7 + 5 - 07); 7 + 5 - a, 

^ + S-/3,l-a-/3 + ^ + 5,5;l-y), (A4) 

The solution in the neighborhood oi y = can be written as a linear combination of last 
two Heun's functions [14] 

H{a, b; a, (3, 7, 5; y) = DiH{l — a,—b — a(3; a, (3,1 + a + (3 — j — S,d;l — y) 

+ D2(i - yy+^-^-^ 

X H{1 -a,-b-al3-(-f + 5-a-l3)(j + 5-aj); 

7 + 5-a,7 + 5-/3,l-a-/3 + 7 + 5,5;l-7/), (A5) 

where the constants are given 

Di = H{a,b;a,l3,-f,5;l) , 

D2 = H{a, b - 07(7 + 5 - a - 7 + 5 - a, 7 + 5 - /3, 7, 5; 1) . (A6) 

The following identity links the arguments {(3, 7, 6) to (/3 + 1, 7 + 1, 5 + 1), respectively, 

{-fap + b)H{a, b - a;a, P + + 1,S + l;y) 
= a-fH{a, b- a, /3, 7, 5- y) + 07(7/ - l)^H{a, 6; a, /3, 7, 5; y) . (A7) 

Finally, in the limit of a — > 00, Heun's function turns into a hypergeometric function [14] 

H{a, aA; a, /3, 7, S; y) =,^00 2^1 Q(a + - 5) + \j \^{a + 13 - 5)f + A , 

1(« + /3 - 5) - y^[l(« + /3-5)]2 + A; 7; y^ . (A8) 

with7 7^-n(n = 0, 1,2,...). 
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FIG. 1: The reflection and transmission coefficients in the case of position-dependent mass 

for mo = 0.1, mi — 10,5 — 5,Vo — 5. 
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(a) reflection and transmission coefficients for mi = 0.01, (5 = 3, Vb = 1 and E = 0.05. 
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(b) reflection and transmission coefficients for mo = 1, (5 = 5, Vb = 1 and E = 0.1 . 



FIG. 2: Reflection and transmission coefficients in tlie case of position-dependent mass. 
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FIG. 3: The unnormalized wave functions in the case of constant mass for 
m = 1, 5 = 2, ^0 = 0.5, E = -m/10 (sohd hne) and for m = 2, 5 = 2, Vb = 0.5, E = -m/10. 
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(a) reflection and transmission coefficients for m = 0.5, S = 5,Vo = 1. 
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(b) reflection and transmission coeflScients forTO=l,5 = 5,Vb = l. 
FIG. 4: Reflection and transmission coefficients in tlie case of constant mass. 
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